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We study the high density region of QCD within an effective model obtained in the frame of the 
hopping parameter expansion and choosing Polyakov-type loops as the main dynamical variables 
representing the fermionic matter. This model still shows the so-called sign problem, a diffi- 
culty peculiar to non-zero chemical potential, but it permits the development of algorithms which 
ensure a good overlap of the simulated Monte Carlo ensemble with the true one. 
We review the main features of the model and present results concerning the dependence of var- 
ious observables on the chemical potential and on the temperature, in particular of the charge 
density and the Polykov loop susceptibility, which may be used to characterize the various phases 
expected at high baryonic density. In this way, we obtain information about the phase struc- 
ture of the model and the corresponding phase transitions and cross over regions, which can be 
considered as hints about the behaviour of non-zero density QCD. 
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1. Introduction 

The exploration of the phase diagram of QCD at non-zero baryon density is a challenging 
and interesting problem. In particular, it has been emphasized that quark matter at extreme high 
density may behave as a color superconductor and it is also expected that the phase diagram in 
the temperature-density plane shows multiple phases separated by various transition lines but very 
little is known about their exact position and nature. Lattice calculations, using different approaches 
that try to evade the sign problem generated by the non-zero chemical potential, have been mostly 
implemented at high temperature and small baryon density, where they agree reasonably well with 
each other. In this region there is good evidence for the presence of a crossover instead of a sharp 
deconfining transition. At large \i (baryon density), however, there are only few numerical results 
which need to be corroborated by using different methods. See [|l|] for a review. 

The purpose of this work is to get further insight about the phase structure of high density, 
strongly interacting matter using an effective model. In the spirit of the \i = quenched ap- 
proximation a 'non-zero density quenched approximation' for ji > based on the double limit 
M — > oo 5 — > oo 5 £ = exp (ji — InM) : fixed f|] has been considered. This implements a static, 
charged background, which influences the gluonic dynamics [||, |J. The present model [||] repre- 
sents a systematic extension of the above considerations: the gluonic vacuum is enriched by the 
effects of dynamical quarks of large (but not infinite) mass, providing a large net baryonic charge. 
In [^[0, §] we explore the phase structure of the model, as a first step in understanding the properties 
of such a background. 

2. QCD at Large Chemical Potential 

We use the QCD grand canonical partition function with Wilson fermions at ji > defined as: 

&{p,K,y G ,yF,H) = J[DU]z- s ^^ u »^ F {K,y F ^,{U}), (2.1) 
S G (p, 7G ,{U}) = -^ReTr (1 £ ^. + 7g £/> 4 ) , (2.2) 
K,n) = J [DU]c- s ^ u ^^ F (K,^ {U}) , & f (k,h,{U}) = DetW(K,n,{U}) , (2.3) 

3 

w ff = d ff [i -KfY, {r+iUiTi+r.iTfu*) - Kf (^fr +4 u 4 T 4 + e -^r_ 4 T:u;)] , 

r ±M = l±y M , ft = i = 1 >* = 2 (M + 3 + cosh A t) = 2(M + 4) ' 

where we used Sg for Wilson's plaquette (P) action and used a certain definition of the Wilson 
term in W. Here M is the 'bare mass', Mq the bare mass at \i = 0, / is the flavor index, 
the link variables and lattice translations. A non-zero physical temperature T is introduced as 
aT = where Yphys is the physical cutoff anisotropy defined by an appropriate renormalization 
of the coupling anisotropics, and 7V T the 'length' of the (periodic) temporal lattice size. 

At large ji the behavior of QCD quantities may however be dominated by certain factors in 
the fermionic determinant which lead to a simpler model that is actually easier to simulate. The 
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Figure 1: Periodic lattice, loops, temporal gauge. In the maximal temporal gauge also the links of the basis 
line are fixed to 1 up to the rightmost one. 



model we study is based on an analytic expansion of QCD (the hopping parameter expansion) up to 
second order in k (see H for details) and its main ingredient are Polyakov-type loops (see fig. [I]), 
capturing the effect of heavy quarks with low mobility. The model still has a sign problem, being 
the fermionic determinant complex at ji ^ 0, but being based on the variables which are especially 
sensitive to the physics of dense baryonic matter it allows for reweighting algorithms which ensure 
a good overlap of the Monte Carlo ensemble with the true one. 

We use the Wilson action and Wilson fermions within a reweighting procedure. The updating 
is performed with a local Boltzmann factor which only leads to a redefinition of the "rest plaquette": 



Plaq x ^< 



The weight (global, vectorizable) is 



(2.4) 



w{{U}) = Y\^ V i-2CReTx [^+K 2 £^ J Wj 2] ({[/}), (2.5) 
x I L itf J J 



such that the 'Boltzmann factor' becomes, 



(2.6) 



wB = B = Yl e^ ReTrPlaq ^ 2 \{U}) . 

Plaq 

Averages are calculated by reweighting according to, B = Bowo and (O) = • 

We have employed the Cabibbo-Marinari heat-bath procedure mixed with over-relaxation. 
This updating already takes into account part of the ji > effects and the generated ensemble 
can thus have a better overlap with the true one than an updating at ji = 0. One can also use an 
improved Bo, to be taken care of by a supplementary Metropolis check. Notice that extracting a 
factor like Bq may also improve convergence of full QCD simulations at ji > 0. 
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Figure 2: Comparison with strong coupling at /3 = 3 (left plot) and j3 = 5 (center plot), 4 4 lattice. Full 
symbols denote ReP, empty symbols ReP*, the lines show the corresponding strong coupling results. On the 
right is reported the mean field phase diagram (abscissa /jl, ordinate y = N T a T). 



We measure several observables under the variation of \i and T, to check the properties of the 
different phases for small T and large jl. In the following we specialize to N c = 3. The observables 
are: the Polyakov loop ((P)) and its susceptibility (Xp) 

G x G x y 

the (dimensionless) baryonic number density rip = E/ yf* with the corresponding susceptibility. 
To check the character of the conjectured third phase we also measure (but we do not report the 
corresponding results here) the spatial and temporal plaquettes, the topological charge, topological 
charge susceptibility and the diquark-diquark correlators (see [||] for details). 

The simulations are mainly done on lattice 6 4 for rif = 1,3 degenerate flavors (any mixture of 
flavors can be implemented). The K dependence has been analyzed in Here we set K = 0.12 
(rather "small" bare mass Mo = 0.167) which drives the 1/M 2 effects in the baryonic density to 
about 50%. The task we have set to ourselves is primarily to explore the phase structure of the 
model at large chemical potential and "small" temperature and we accordingly vary ji and j8. 



3. Results and Discussions 



As a first orientation about the behavior of the model we consider the strong coupling/hopping 
parameter expansion, which also serve as a check of the Monte Carlo results. In fig. || we compare 
the Monte Carlo results of the Polyakov loop and its adjoint on 4 4 and 6 4 lattices, for K = 0.12, 
one flavor and different values of j8, with and P*^ (where [2] means second order in the strong 
coupling expansion, see ^ for details). The agreement is good for the 4 4 lattice and j8 = 3, while 
for j8 = 5 there are already significant deviations which show strong effects at large ji even at 
moderate j8 that may indicate possible phase transitions. But the agreement between Monte Carlo 
and strong coupling results is sufficient to validate the simulations. 

We also performed mean field calculations using a temporal gauge fixing [§] and introducing 
two different mean field variables for the spatial component (u) and the temporal component (v) of 
the gauge field. They give some qualitative insight into the phase structure of the model to which 
Monte Carlo simulations can be compared. In fig. [2| we give an illustrative example, taken with 
j8 = 4 and 7V T = 6. It shows a large 'confinement' region for small T and ji corresponding to the 
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Figure 3: Left: Landscape of the baryonic density. Right: Landscape of the baryon density susceptibility. 
The color scale (right) is based on log 10 (ft#). 

trivial fixed point mentioned above with both mean fields u and v vanishing. For larger T or ji one 
crosses into a deconfined regime with both mean fields w,v > 0. In the lower right corner there 
appears in addition an intermediate phase with u = 0, v > 0. The field v is close to its maximal 
value 1 wherever it is nonzero, whereas u has smaller, varying values, depending on the region. 

The algorithm works reasonably well over a large range of parameters even at small temper- 
ature. The model permits to vary \i, /c, j8 as independent parameters and it is reasonably cheap 
to measure various correlations. In fig. || we show the "landscape" of the real part of the baryon 
density hb (while the imaginary part is compatible with zero inside the statistical errors). The main 
variation is an exponential growth with ji indicating that we do not see yet saturation effects. This 
masks to a certain extent the finer structure. A clearer view of the situation is provided by looking 
at the "landscape" of the susceptibility of the baryon density (fig |3j) . A ridge is clearly visible, 
highlighted by a dashed black line. A second line (dashed) is explained later. We found it therefore 
advantageous to look at the Polyakov loops and their susceptibility. In fig. || we show the Polyakov 
loop susceptibility vs. j8 at fixed ji and on the bottom vs. fi at fixed j8 and, in fig || (right) the 
corresponding landscape. The plots of the Polyakov loop susceptibility show quite clearly maxima 
indicating possible transitions or crossovers. In the landscape fig. ||, one of these maxima shows 
up as a well defined ridge, indicated by a dashed black line. It shows only a moderate slop in /i, 
which explains why the maxima are more pronounced when we vary j8 at fixed \i than vice versa. 

The broadening of this ridge at small \i as well as of the maximum in Fig. |] is responsible 
for the loss of a sharp transition signal at small ji. These figures clearly show that the transition at 
fixed ji = 0.50 is less steep than the one at ji = 0.80. Presumably at ji <~ 0.6 we are dealing with 
a crossover, whereas at large ji the signal is more compatible with a real phase transition. Notice 
that changing j8 at fixed /I, we cross the transition line at a more oblique angle at smaller /I, but 
the broadening of the ridge and loss of a transition signal is a genuine effect, as can be seen from 
fig. ||. A second ridge branching off from this main ridge at large /I, highlighted by a dotted line is 
suggested by looking at the level lines in fig. |5| and corresponds to the second maximum suggested 
at large \i in fig. |J This may indicate the appearance of the new phase at large ji and small T/T c 
discussed above. 

We used the results for the Polyakov loop susceptibility to estimate the possible position of the 
transition points in the j8 vs. ji plane (see ^ for details); to go half way toward a possible physical 
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Figure 4: Polyakov loop susceptibility vs. j3 at fixed ji (top) and vs. ji at fixed j3 (bottom). 




|X phys /T c = 6|iT/T c 



Figure 5: Left: the phase diagram in the j3 (or T/T c ) - ji p hys/T c QCD plane. The dotted straight lines 
correspond to constant /I, the dashed ones to constant j3 . Right: landscape of the Polyakov loop susceptibility 
(top) and its 3d view (bottom). The color scale is based on log 10 (^p). 
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interpretation the positions determined in this way are indicated by the blobs in the diagram T/T c 
vs. Uphys/Tc of fig. ||, where fl p hy S = l^/ a (P) = N T jiT. The shaded blobs correspond to the 
rather unambiguous 'deconfining' signal observed for ji >~ 0.6 (j8 <~ 5.72). The 'transition' line 
suggested by this signal starts at the lower point A on the figure, located at j8 ~ 5.55, \i ~ 0.88, i.e., 
with our rough estimation Hphys/T c — 2.4, T/T c ~ 0.45 (below which we could no longer obtain 
reliable data) and ends at the point B located near j8 ~ 5.72, ji ~ 0.6, i.e., with our rough estimation 
Uphys/Tc — 2.3, T/T c ~ 0.65. Above this point the signal becomes ambiguous. 

The picture emerging from the data is thus the following: for ji < 0.5 — 0.6 (jiphys/T ~ 3) 
there is only a broad crossover, while for 0.6 < ji < 0.9 (3.6 < li p h ys /T < 5.3) there is evidence of 
a sharper crossover or transition at a value ji c depending on j8. Moreover, for ji ~ 0.9 there is some 
evidence of the presence of the second transition even though this evidence is much weaker than 
the other one because at larger values of ji the fermion determinant strongly oscillates and, indeed, 
the usual sign problem manifest its effects. 

To summarize our results, the phase structure found by the numerical simulations for rif = 3 
is shown in fig. ||. The signal for the deconfining transition (or narrow crossover) on the line 
connecting A and B is rather good and it also appears that at small ji (above B) the transition is 
smoothed out in accordance with the expectations from full QCD simulations [p], fj. A second 
transition at large ji could only be identified tentatively. In this region, the diquark susceptibility 
grows strongly This region needs further study to reach a conclusion, but it is interesting that 
the general picture shows qualitative agreement with the one found in the mean field approximation. 

We can consider this model as an evolved 'quenched approximation' in the presence of charged 
matter. Then this study would give us information about the modified gluon dynamics of the SU(3) 
theory in this situation. It would then be natural to think of it as providing a heavy, dense, charged 
background for propagation of light quarks and calculate light hadron spectra and other hadronic 
properties under such conditions. This could also help fixing a scale controlling the behavior of the 
light matter. Work in progress goes in this direction. 
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